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Measurements  have  been  made  of  first-motion  amplitude  and  of  the  period 
and  amplitude  of  the  second  peak-to-trough  excursion  (termed  the  £ phase)  for 
short-period  vertical  P-wave  signals  from  approximately  60  Nevada  Test  Site 
(NTS)  explosions.  They  were  recorded  at  the  Long-Range  Seismic  Measurement 
(LSRM)  and  Seismic  Data  Collection  System  (SDCS)  stations  in  Houlton,  Maine 
(HNME) , and  Red  Lake,  Ontario  (RKON) . Almost  all  events  selected  were  deto- 
nated in  saturated  tuff  or  rhyolite.  Care  was  taken  to  avoid  gross  errors 
and  to  minimize  effects  of  apparent  minor  perturbations  of  the  waveforms,  - — 
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The  £-phase  data  are  plotted  as  a function  oft  Ms  and  are  interpreted  to 
show  amplitude  variations  of  several  tenths  in  magnitude  with  respect  to  an  Ms; 
yield  curve  with  slope  of  1.0.  These  amplitude  variations  are  due  independently 
to:  (1)  the  influence  of  yield  tin  corner  frequency;  (2)  interference  with  the 

pi’  phase;  and  (J)  geographic  location  within  NTS.  Several  tenths  of  a second 
variation  in  period  can  also  he  traced  to  the  above  3 causes.  The  variation 
with  geographic  location  may  be  due  to  variation  of  the  upper  2 km  of  structure 
at  NTS  or  to  multipathing  and  ray  divergency  because  of  more  deeply  seated 
structures. 

V 

in  general  the  data  supports  the  conclusion  that  the  nq,:  yield  curve  begins 
to  bend  significantly  at  approximately  200  kt,  as  predicted  by  the  granite 
reduced  displacement  potential  (Kl)l’)  of  von  Seggern  and  Hlandford.  The  curva- 
ture may  not  have  been  detected  earlier  at  NTS  due  to  scatter  in  the  data  from 
the  three  effects  mentioned  above.  Another  possibility  shown  in  tills  report  is 
that  near  the  1000-kt  region  the  slope  of  log  (c/T)  versus  log  (yield)  is  closer 
to  1.0  than  is  the  slope  of  log  (O  versus  log  (yield). 

The  scatter  of  the  data  is  greater  at  KKON  than  at  HNME,  presumably  re- 
flecting the  greater  variability  of  signals  at  regional  distances  and  the 
shallow  takeoff  angles.  Network  averages  of  teleseismic  stations  would  certainly 
show  less  scatter.  Despite  the  greater  amount  of  scatter  in  single-station  data, 
we  feel  that  the  general  trends  are  clear  enough  to  support  the  conclusions 
reached.  Use  of  single-station  data  is  useful  because  there  are  fewer  "hidden" 
effects  of  data  selection. 
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ABSTRACT 


Measurements  have  been  made  of  first-motion  amplitude  and  of  the  period 
and  amplitude  of  the  second  peak-to-trough  excursion  (termed  the  £-phase)  for 
short-period  vertical  P-wave  signals  from  approximately  60  Nevada  Test  Site 
(NTS)  explosions.  They  were  recorded  at  the  Long-Range  Seismic  Measurement 
(LSRM)  and  Seismic  Data  Collection  System  (SDCS)  stations  in  Houlton,  Maine 
(HNME),  and  Red  Lake,  Ontario  (RKON).  Almost  all  events  selected  were  deto- 
nated in  saturated  tuff  or  rhyolite.  Care  was  taken  to  avoid  gross  errors 
and  to  minimize  effects  of  apparent  minor  perturbations  of  the  waveforms, 
which  cause  zero-crossings,  the  so-called  "false  cycles". 

The  jC-phase  data  are  plotted  as  a function  of  Ms  and  are  interpreted  to 
show  amplitude  variations  of  several  tenths  in  magnitude  with  respect  to  an 
Ms:  yield  curve  with  slope  of  1.0.  These  amplitude  variations  are  due  inde- 
pendently to:  (1)  the  influence  of  yield  on  corner  frequency;  (2)  interfer- 
ence with  the  pP  phase;  and  (3)  geographic  location  within  NTS.  Several 
tenths  of  a second  variation  in  period  can  also  be  traced  to  the  above  3 
causes.  The  variation  with  geographic  location  may  be  due  to  variation  of 
the  upper  2 km  of  structure  at  NTS  or  to  multipathing  and  ray  divergency 
because  of  more  deeply  seated  structures. 

In  general  the  data  supports  the  conclusion  that  the  m^:  yield  curve 

2 

begins  to  bend  significantly  at  approximately  200  kt,  as  predicted  by  the  w 
granite  reduced  displacement  potential  (RDP)  of  von  Seggern  and  Blandford. 

The  curvature  may  not  have  been  detected  earlier  at  NTS  due  to  scatter  in  the 
data  from  the  three  effects  mentioned  above.  Another  possibility  shown  in 
this  report  is  that  near  the  1000-kt  region  the  slope  of  log  (c/T)  versus  log 
(yield)  is  closer  to  1.0  than  is  the  slope  of  log  (c)  versus  log  (yield). 

The  scatter  of  the  data  is  greater  at  RKON  than  at  HNME,  presumably  re- 
flecting the  greater  variability  of  signals  at  regional  distances  and  the 
shallow  takeoff  angles.  Network  averages  of  teleseismic  stations  would  cer- 
tainly show  less  scatter.  Despite  the  greater  amount  of  scatter  in  single- 
station  data,  we  feel  that  the  general  trends  are  clear  enough  to  support  the 
conclusions  reached.  Use  of  single-station  data  is  useful  because  there  are 
fewer  "hidden"  effects  of  data  selection. 
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Figure  Title  Page 

1 Locations  of  events  discussed  in  this  report  at  the  No-  lb 

vada  Test  Site.  F.vents  are  indicated  by  numbers,  upper 

and  lower  case  letters,  and  other  common  symbols.  The 
event  associated  with  each  symbol  is  given  in  Table  1. 

Symbols  for  source  regions  are  indicated,  as  well  as 
azimuth  to  RKON.  The  outline  of  the  Silent  Canyon  Cal- 
dera, as  found  by  Orkild  et  al.  (19b9),  is  drawn  around 
the  Pahute  Mesa  events;  and  the  surface  trace  of  Yucca 
Fault,  as  mapped  by  Fernald  et  al.  (19b9),  is  drawn 
through  the  Yucca  Flats  events.  Height  above  the  water 
table  is  indicated  in  meters  for  the  non-RAINlER  events 
above  the  water  table. 

2 Waveforms  of  the  events  I'.AMEMBERT  (z)  and  KASSERI  (0)  as  18 

seen  at  HNME  and  RKON.  Measurements  a,  b,  and  c are  in- 
dicated. For  reasons  discussed  in  the  text,  the  defini- 
tion for  £ is  one-half  cycle  later  than  that  usually 

used.  Note  that  small,  high-frequency  "false  cycles" 
are  ignored  in  measurements  of  both  amplitude  and  period. 

Ja  Theoretical  amplitude-yield  curves  for  t*  ■ 0,  0.1,  0.2,  19 

0.4,  O.b;  granite,  amplitude  of  first  motion,  with  no 
surface  reflection.  From  Rlandford  (197b). 

Jb  Theoretical  amplitude-yield  curve  for  t*  ■ 0,  0.1,  0.2,  20 

0.4,  O.b;  granite,  amplitude  of  first  motion,  with  sur- 
face reflection.  pP  delay  equal  to  0.12Y1'  * sec,  with  Y 
in  kt . From  Blandford  (197b). 

3c  Theoretical  amplitude-yield  curves  for  t*  ■ 0,  0.1,  0.2,  21 

0.4,  O.b;  granite,  one-half  maximum  peak-to-t rough  am- 
plitude of  signal,  with  no  surface  reflection.  From 
Blandford  (197b). 

Id  Theoretical  amplitude-yield  curves  for  t*  ■ 0,  0.1,  0.2,  22 

0.4,  O.b;  granite,  one-half  maximum  peak-to-t rough  am- 
plitude of  signal,  with  surface  reflection.  pP  delay 
equal  to  0.12Y1'  ' second,  with  Y in  kilotons.  From 
Blandford  (197b). 

3e  Theoretical  amplitude-yield  curves  for  t*  * 0,  0.1,  0.2,  23 

0.4,  O.b;  granite,  one-half  maximum  peak-to-trough  am- 
plitude of  signal,  corrected  for  instrument  response  at 
measured  period,  T,  and  divided  by  T;  with  no  surface  re- 
flection. From  Blandford  (197b). 
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it  Theoretical  ampl ltude-y leld  curves  for  t * - 0,  0.1,  0.2,  24 

0.4,  0.t>;  granite,  one-halt  maximum  peak-to-t  rough  am- 
plitude tor  signal,  corrected  toi  instrument  response  at 
measured  period,  T,  and  divided  hy  T,  with  surtace  retleet- 
tion.  pi'  delay  equal  to  0.12Y1/I  sec,  with  Y in  kt . 

From  Hlandtord  (ld/b), 

4a  RKUN  Station:  l.ogfa)  as  a function  ot  Ms.  The  dashed  27 

line  Is  the  least -squares  line  ot  slope  ^.0  through  all 
the  data.  The  solid  lines  are  t rom  Figure  la  with  t*  ■ 

0.4,  with  each  line  vertically  adjusted  hy  eve  so  t hev 
respectively  best  tit  the  west  1’ahute  data  and  all  other 
dat  a. 

4b  Residuals  with  respect  to  the  dashed  line  In  Figure  4a  as  jg 

a t unction  ot  latitude  and  longitude  at  NTS.  Tiro  residuals 
at  West  1’ahute  have  been  contoured.  F.vents  may  be  identi- 
t led  by  comparison  with  Figure  l and  Table  l.  Solid  dots 
are  events  lacking  any  data. 

la  RKUN  Station:  l.og  (ji)  as  a function  ot  logic).  The  12 

dashed  line  is  the  least-squares  line  ot  slope  1.0 
through  all  the  data. 

5b  RKUN  Station:  Residuals  ot  log(a)  with  respect  to  the 

dashed  line  in  Figure  5a  as  a t unction  ot  latitude  and 
longitude  at  NTS.  Kvents  are  identified  in  Table  l and 
Figure  1. 

5c  KK0N  Station:  Residuals  ot  log(a)  with  respect  to  the  14 

dashed  line  lit  Figure  5a  as  a (unction  ot  scaled  2-wav 
pi’  travel  time.  The  solid  line  is  hand-drawn  as  a sub- 
jective estimate  ot  the  main  data  trends. 

5d  RK0N  Station:  Wavetorms  ot  the  events  SCOTCH,  COMMOOORK , IS 

RF.X,  and  TYRO  to  show  the  variety  ot  possible  signals. 

ba  ^ as  a function  of  M at  KKON.  The  dashed  line  is  the  least-  17 

squares  lines  ot  slope  1,0  through  all  the  data.  The  solid 
curved  lines  are  t rom  Figure  1c  with  t*  ■ 0.4  and  with  the 
vertical  ordinate  adjusted  by  eye  separately  to  tit  tire  West 
1’ahute  and  all  other  data.  Note  that  the  Hast  I’ahute  data  also 
fit  this  latter  curve  (in  contrast  to  the  data  in  Figure  4a). 

bb  Residuals  ot  log  (e^)  with  respect  to  the  dashed  line  in  18 

Figure  ba  as  a function  ot  latitude  and  longitude  at  NTS. 

Kvents  can  be  Identified  bv  comparison  with  Figure  1 and 
Table  l. 
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KKON  Station:  l-oglc),  with  c corrected  tor  pi’,  as  a id 

t’ unction  ot  Ms.  pi’  correct  Ion  is  given  hv  the  solivl 

line  in  Figure  Sc.  I'he  dashed  line  is  the  least- 

squares  line  ot  slop*'  1.0  through  all  t lie  data.  The 

solid  curved  lines  are  from  Figure  Jc  with  t * * 0.4  and 

with  .he  vertical  ordinate  adjusted  by  eye  separately  to 

tit  the  West  i'ahute  and  all  other  data. 

KKON  Station:  Period  ot  c as  a function  ot  Ms.  The  40 

dashed  1 in*'  is  the  1 east -squares  regression  on  M . . 

KKON  Station:  Kesiduals,  in  tenths  ot  a seconds,  *>t  t ho  42 

c period  with  respect  to  the  dashed  line  in  Figure  8a  as 
a function  of  latitude  and  longitude  at  NTS.  "A"  cor- 
responds to  0.  I to  0. 2 set'  below  the  lino;  "l"  corres- 
ponds to  0.  1 to  l)..’  sec  above  the  1 in*'.  Fvents  are 
identified  in  Table  l and  Figure  l. 

KKON  Station:  Residuals,  in  tenths  ot  a second,  ot  the  s' 

c period  with  respect  to  the  dashed  line  in  Figure  8a  as 
a function  ol  scaled  2-way  pi’  travel  time.  I’he  solid  line 
has  been  subjectively  drawn  to  point  out  the  trend  in  the 
data.  The  RKON  waveforms  tor  2 events  at  either  end  ot 

the  trend,  PILEDRIVKR  ())  and  Q)BB1.KK  (P)  are  Included  in 

the  upper  left-hand  corner. 

c/1'  with  £ corrected  for  Instrument  response  at  period  1 as  ^ 
a function  of  M at  KKON.  The  dashed  line  is  the  least- 
squares  line  of  slope  1.0  through  all  the  data.  The 
solid  curved  lines  are  from  Figure  le  with  t*  - 0.4  and 
with  the  vertical  ordinate  adjusted  bv  eve  separately  to 
fit  the  west  I’ahute  and  all  other  data. 

Kesiduals  with  respect  to  the  dashed  l in*-  in  Figure  8.*  as  4o 

a function  ot  latitude  and  longitude  at  NTS.  Fvents  may 
be  identified  bv  comparison  with  Figure  1 and  Table  l. 

a as  a function  ot  M at  HNMF.  The  dashed  line  is  the  47 

least-squares  line  ot  slope  1.0  through  all  the  data. 

The  solid  curved  line  is  from  Figure  3a  with  t*  * O.b  and 
with  the  vertical  ordinate  adjusted  bv  eye  t o tit  t he  data. 
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Title  Pane 

First -mot  Ion  amplitude  residuals  at  HNMK  with  respect  to  4b 

the  dashed  line  in  Figure  10a  as  a function  ot  latitude  and 
longitude  at  NTS.  Events  may  be  identified  by  comparison 
with  Figure  l and  Table  1. 


a as  a function  of  e at  HNMK.  The  dashed  line  is  the  least- 
squares  line  ol  slope  1.0  through  all  the  data. 

Residuals  of  log  (a)  at  HNMK  with  respect  to  the  dashed  SO 

line  in  Figure  11a  as  a function  of  latitude  and  longi- 
tude. 


Residuals  ot  log  (a)  at  HNMK  with  respect  to  the  dashed  line  51 
in  Figure  11a  as  a function  ot  scaled  2-way  pi’  travel  time. 

The  solid  line  is  hand-drawn  as  a subjective  estimate  ol 
the  mean  data  trends. 

£ as  a function  of  M at  HNMK.  The  dashed  line  is  the  32 

least-squares  line  of  slope  1.0  through  all  the  data. 

The  solid  curved  lines  are  from  Figure  3c  with  t*  - O.b 
and  with  the  vertical  ordinate  adjusted  separately  by  eye 
to  fit  the  West  Pahute  and  all  other  data. 

c residuals  with  respect  to  the  dashed  line  in  Figure  12a  51 

at  HNME  as  a function  of  latitude  and  longitude  at  NTS. 

Events  may  be  identified  by  comparison  with  Figure  1 and 
Table  1. 

£ plus  pP  correction  for  £ at  HNMK  as  it  function  of  Ms.  55 

pP  correction  given  by  the  solid  line  seen  in  Figure  lie. 

The  dashed  line  is  the  least-squares  line  of  slope  1.0 
through  all  the  data.  The  solid  curved  lines  are  from 
Figure  3a  with  t*  = 0.6  and  with  the  vertical  ordinate 
adjusted  separately  by  eye  to  fit  the  west  Pahute  and  all 
other  data. 

Period  of  c at  HNMK  as  a function  of  M . The  dashed  line  5h 
is  the  least-squares  regression  on  M^.‘s 

Residuals  of  the  £ period  at  HNMK  in  tenths  of  a second  57 

with  respect  to  the  dashed  line  in  Figure  14a  as  a function 
of  latitude  and  longitude.  "A"  corresponds  to  0. 1-0.2 
seconds  above  the  line.  Events  may  be  identified  by  com- 
parison with  Figure  l and  Table  I. 
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14c  HNMK  Station:  Residuals,  in  tenths  of  a second,  ot  the 

£ period  with  respect  to  the  dashed  line  in  Figure  14a 
as  a function  of  scaled  --way  pP  travel  time.  The  solid 
line  is  hand-drawn  as  a subjective  estimate  ot  the  main 
data  trends. 

ISa  HNMK  Station:  Log^c/T),  with  c corrected  tor  instrument 

response  at  period  V,  as  a function  ot  Ms.  The  dashed 
line  is  the  least-squares  line  of  slope  1.0  through  all 
the  data.  The  solid  line  is  from  Figure  Je  with  t*  » 
0.t>,  with  the  line  vertically  adjusted  bv  eye  so  it  best 
fits  the  data. 

1 Sb  HNMK  Station:  Residuals  ot  log(c/T)  with  respect  to  the 

dashed  line  in  Figure  lba  as  a function  of  latitude  and 
longitude  at  NTS.  Kvents  are  identified  in  Table  1 and 
Figure  l. 
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INTRODUCTION 

l'lu'  pur  post'  of  this  report  is  to  examine  In  detail  the  variation  oi  sin- 
gle-station, short-period  i’  parameters  derived  t rora  measurements  ot  the  Neva- 
da lest  Site  (NTS)  explosion  signals  at  single  stations  (KkON  and  IINMT)  as  a 
function  ot  Ms.  Because  Ms  is  a well-determined  function  ot  vield  (von  Seg- 
gern, 117b),  this  study  is  equivalent  to  an  examination  ot  the  variation  as  a 
function  ot  vield.  Effects  at  single  stations  may  more  clearly  show  el  toots 
of  epicentral  geology.  As  ot  this  writing  (October  1177),  final  computet 
outputs  are  combining  t he  results  in  this  paper  with  those  from  l other  sta- 
t Ions  to  produce  network-average  results.  Preliminary  results  indicate 
little  change  in  the  qualitative  conclusions. 

Several  other  authors  have  already  approached  this  problem.  Among  them 

are  Marshall  et  a l . (I1*/’ 11  and  von  Seggern  (117b),  who  have  discussed  data 

which  confirm  simple  cube-root  scaling  predictions  ot  an  M . : log (y ie Id! 

slope  ot  1.0.  von  Seggern  (117b)  found  that  M from  t hi'  M'S  explosions 

s 

below  the  water  table  .it  Pahute  Mesa  were  0.1  and  0.2  magnitude  units  higher 
than  M values  for  equivalent-'1  ield  shots  below  the  water  table  elsewhere  at 
NTS.  von  Seggern  also  pointed  out  that  the  scatter  of  the  M tyield  data 
points  within  regions  at  NTS  is  very  small  for  events  below  the  water  table. 
He  also  showed,  however,  that  it  onlv  pure  teleseismic  IK  measurements  in  the 
range  of  17  to  21  sec  are  used,  then  for  corresponding  yields  (BOXCAR  and 
Mil. ROW),  M is  0.8  units  larger  for  explosions  at  NTS  than  at  Amchltka. 

Since  from  each  region  these  events  tall  on  a very  tight  M tvield  curve  ot 
slope  1.0,  this  discrepancy  would  apply  to  any  pair  ot  events,  one  at  each 
test  site. 

The  above  findings  show  the  importance  of  m^tyield  relat tonshlps  for  a 

threshold  test-ban  treaty.  Despite  controversy  ovot  the  existence  ot  m^-mag- 

nttude  bias  at  different  test  sites  (Her,  117b),  m^  bias  is  generally  thought 

von  Seggern,  D,  H.  (117b).  M,  versus  yield  of  underground  nuclear  explosions 
at  the  Nevada  Test  Site;  SDAC-TR-76— l l , Teledvne  Coot  cell,  Alexandria, 
Virginia. 

Marshall,  P.  D. , A.  Douglas,  and  J.  A.  Hudson  ( l v> 7 1 ) . Surface  waves  from 
underground  nuclear  explosions;  Nature,  2 1-),  8-1. 

Her,  Z.  A.  (117b).  On  the  existence,  magnitude  and  causes  of  broad  regional 
variations  in  body-wave  amplitude  (magnitude  bias);  Sl)AC-l’R-7b  8,  Tele- 
dvne  t’.eotech,  Alexandria,  Virginia. 
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to  be  less  than  0.5  units  (0.1  < bias  < 0.25),  and  most  workers  accepting  the 
existence  of  bias  also  think  it  can  be  removed  by  careful  attention  to  de- 
tailed corrections  applied  to  the  station  magnitudes  comprising  network  mag- 
nitude averages.  Thus,  m^  may  be  a more  reliable  estimate  of  yield  than  Mg, 
making  m^:yield  relations  worthy  of  more  serious  study. 

Springer  and  Hannon  (1973)  presented  m^ryield  data  for  explosions  at 
Pahute  Mesa.  They  concluded  that  the  m^: log(yield)  slope  is  approximately 
1.0  for  teleseismic  distances,  and  approximately  0.6  for  regional  distances. 
They  attributed  the  difference  to  the  total  effect  of  absorption  which  in- 
creases with  distance,  and  then  they  developed  theoretical  signal  spectra 
which  implied  waveforms  that  would  produce  such  a change  of  slope  with  dis- 
tance. 

Blandford  (1976)  computed  observed  regional  spectral  ratios  for  Pahute 
Mesa  events  ranging  in  magnitude  from  0.7  to  1200  kilotons.  Using  the  RAINIER 
reduced  displacement  potential  (RDP)  for  tuff,  he  found  that  the  spectral  ra- 
tios and  relative  amplitudes  agreed  with  the  predictions  of  theory  outlined 
by  von  Seggern  and  Blandford  (1972).  The  same  conclusion  was  reached  by 
Noponen  (1975)  in  comparing  granite  RDP  to  teleseismic  observations  of  USSR 
and  NTS  explosions  and  by  Lyuke  et  al.  (1976)  for  regional  observations  of 
USSR  explosions.  Blandford,  like  Rodean  (1972),  calculated  theoretical 

Springer,  D.  L. , and  W.  J.  Hannon  (1973).  Amplitude-yield  scaling  for  under- 
ground nuclear  explosions;  Bull.  Seism.  Soc.  Am. , 63,  477-500. 

Blandford,  R.  R.  (1976).  Experimental  determination  of  scaling  laws  for  con- 
tained and  cratering  explosions;  SDAC-TR-76-3,  Teledyne  Geotech,  Alexan- 
dria, Virginia.  DDC:  ADA  03635. 

von  Seggern,  D.  H.,  and  R.  R.  Blandford  (1972).  Source  time  functions  and 
spectra  for  underground  nuclear  explosions;  Geophys.  J_.  II.  Astr.  Soc. , 

31,  83-97. 

Noponen,  I.  (1975).  Compressional  wave  power  spectrum  from  seismic  sources; 
Institute  of  Seismology,  University  of  Helsinki,  ISBN  951—45-0538-7, 
Contract  AFOSR-72-2377  Final  Report. 

Lyuke,  E.  I.,  S.  K.  Daragan,  and  V.  E.  Peregontseva  (1976).  Forecasting  the 
seismic  wave  spectra  of  large  underground  detonations  from  the  spectra 
of  small  preliminary  detonations;  Phys.  Solid  Earth  (from  Izvestiya), 
12(2),  103-109. 

Rodean,  H.  C.  (1972).  Nuclear-Explosion  Seismology;  U.S.  Atomic  Energy  Com- 
mission, Division  of  Technical  Information,  Oak  Ridge,  Tennessee. 
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waveforms  and  the  associated  m^tyield  curves.  Using  the  tuff  potential, 
Blandford  found  that  absorption  alone  could  not  explain  the  near-unlty  m^: 
yield  teleseismlc  slope.  He  concluded  that  two  RDP's  calculated  from 
measurements  at  angles  of  departure  appropriate  to  regional  and  teleseismlc 
distances  could  not  be  the  same. 

This  report  shows  that  the  granite  RHP  does  fit  the  NTS  teleseismlc  data, 
von  Seggern  and  Blandford  (1972)  and  Blandford  (1976),  using  first-motion 
data  to  avoid  pP  Interference  effects,  detected  the  decrease  In  slope  implied 
by  the  granite  RHP  for  events  LONGSHOT,  MILROW,  and  CANNIKIN,  detonated  at 
Amchitka. 

Mueller  and  Murphy  (1971)  obtained  emprical  scaling  results  for  events 
as  a function  of  depth  that  Implied  a fairly  strong  departure  from  pure  cube- 
root  scaling  (variation  of  the  RHP  as  a function  of  yield) as  depth  of  the 
event  changed.  Note,  however,  that  material  properties  of  the  medium  varied 
with  depth  for  many  of  the  events  used  in  their  regression  analyses  to  de- 
termine empirical  factors.  Using  this  theory,  Murphy  (1976)  presented  theo- 
retical spectral  ratios  for  the  events  RKX  and  SCOTCH,  which  agree  witli  his 
data.  However,  regional  spectral  data  for  these  events  presented  by  Bland- 
ford  (1976)  agreed  with  pure  cube-root  scaling,  and  disagreed  with  Murphy's 
data.  A possible  hypothesis  for  resolving  the  conflict  in  the  experimental 
data  is  to  assume  large  relative  values  of  low-frequency  earth  noise  for 
Murphy's  low-magnitude  RKX  event  data,  and  equally  high  millimicron  noise, 
due  to  high-frequency  system  noise,  for  his  large-magnitude  SCOTCH  event. 
Blandford  (1976)  recognized  and  discussed  such  data-biasing  effects.  The  ef- 
fect of  these  noise  biases  on  the  SCOTCH/RKX  spectral  ratio  is  to  raise  the 
ratio  at  low  frequencies,  and  lower  it  at  high  frequencies.  This  difference 
is  exactly  that  observed  between  the  data  presented  by  Blandford  (1976)  and 
the  data  of  Murphy  (1976).  The  theories  of  the  two  authors  differ  in  the 
same  way  as  do  their  data. 

This  report  treats  a greater  volume  of  NTS  data  than  was  available  to 
Springer  and  Hannon  (1973)  because  it  utilizes  the  Special  Data  Collection 
System  (SOCS)  data  from  stations  RKON  and  HNMK.  Care  is  taken  to 

Mueller,  R.  A.,  and  J.  R.  Murphy  (1971).  Seismic  characteristics  of  under- 
ground nuclear  detonations:  Part  I,  Seismic  spectrum  scaling;  Bull. 

Seism.  Soc.  Am. , 61,  1675-1692. 
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measure  at  each  station  a common  cycle  for  each  event.  This  is  accomplished 
by  tracing  each  signal  and  picking  the  cycle  while  all  the  data  are  displayed 
for  comparison.  Many  potential  gross  data  errors,  seemingly  inevitable  in 
projects  of  this  size,  were  removed  by  careful  reexamination  of  the  original 
waveforms  and  calibration  logs  when  plots  of  the  data  revealed  outlying 
points.  The  data  have  been  stored  in  a computer  file  and  plots  are  produced 
by  computer. 


EVENT  SELECTION  AND  MEASUREMENTS 

Events  selected  for  study  are  listed  In  Table  I.  In  general  they 
consist  of  explosions  with  yields  greater  than  about  10  kt  (HANDCAR,  12kt, 
and  HARDHAT,  5.9  kt , are  among  the  smaller  events)  and  detonated  below  the 
water  table.  Several  events  from  above  the  water  table  are  included,  but  are 
within  a 1-cavity  radius  of  the  water  table.  Also  Included  are  events  ill 
Ranler  Mesa  that  are  above  the  water  table,  but  are  thought  to  be  In  or  near 
layers  of  perched  water  (Kamspott  and  Howard,  1975).  A few  events  well  above 
the  water  table  are  Included  because:  (1)  either  they  are  events  with  shot 
reports  prepared  using  the  SDCS  (EDAM,  MARCH);  (2)  they  are  in  dolomite  (HAND- 
CAR,  KANKAKEE);  or  (3)  they  are  of  special  Interest  for  other  reasons  (DURYEA, 
M1NC  VASE) . . 

Most  of  the  data  in  Table  1 are  taken  from  Springer  and  Kinnaman  (1971). 
Data  for  events  after  1970  have  been  taken  from  data  synopsis  sheets  that  the 
ll.S.  Energy  Research  anil  Development  Agency  (ERDA)  issued  for  each  of  the  re- 
cent tests.  For  several  of  these  reports,  the  average  sound  velocity  from 
shot-point  to  surface  was  estimated  by  vibroseis  or  similar  technique.  it  Is 
these  values  which  are  reported  in  Table  I.  For  events  before  1970,  the 
average  values  for  each  region  (Pahute  Mesa,  2.9  km/sec;  Ranler  Mesa,  2.6  km/ 
sec;  Yucca  Flats,  1.7  km/sec)  reported  by  Ramspott  and  Howard  (1975)  were 
used.  Exceptions  were  made  in  the  case  of  dolomite  and  granite,  which  used 
velocities  of  2.7  and  5.5  km/sec,  respectively.  The  average  regional 
values  are  approximations,  and  though  they  probably  should  not  be  used 
for  detailed  comparison  of  individual  predictions  with  observations, 
they  may  be  adequate  for  the  predominantly  statistical  conclusions  reached  in 
this  study.  The  symbol  for  each  event  (st*e  Table  I)  is  the  same  as  plotted 
on  the  map  (Figure  1),  showing  the  relative  positions  of  the  events.  Small 
perturbations  have  been  made  in  the  locations  of  a few  events  at  Yucca  Flats  so 
Lueir  symbols  would  not  overlap,  but  their  relative  positions  have  been 

Pamspott,  L.  D. , and  N.  W.  Hovard  (1975).  Average  properties  of  nuclear  test 
areas  and  media  at  the  USERDA  Nevada  Test  Site;  UCRL-51948,  Lawrence 
Livermore  Laboratory,  Livermore,  California. 

Springer,  1).  L.  and  R.  1..  Kinnaman  (1971).  Seismic  source  summary  for  ll.S. 
underground  nuclear  explosions,  1961-1970;  Bull.  Seism.  Soc . Am. , 61 , 
1073-1098. 
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TABLE  I (Continued) 


SURFACE  BOUNDARY  OF 
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maintained.  For  example,  if  event  A is  east  of  event  B on  the  ground,  it  is 
also  east  of  event  B on  the  map.  Figure  1 is  also  included  as  a transparent 
overlay  in  the  envelope  attached  at  the  back  of  this  report. 
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Figure  2 gives  the  definitions  of  the  measurements  _a,  jj,  and  £,  made  for 
this  study.  Note  that  the  definition  for  measurement  £ is  one  swing  later 
than  what  is  used  in  most  other  studies.  This  definition  was  chosen  so  it 

more  fully  reflected  the  influence  of  pP.  In  almost  every  case,  this  £ 
measurement  is  also  equal  to  the  maximum  peak- to-t rough  measurement  in  the 
first  4 cycles  used  for  calculating  m^.  Note  in  Figure  2 how  "false  cycles" 
are  ignored  in  making  the  measurements  of  £ and  T at  HNME  and  RKON  for  event 
CAMEMBERT.  The  period  T for  the  £ measurements  is  defined  as  the  interval 
between  zero-crossings  on  either  side.  Note  that  at  HNME,  the  £ measurement 
gives  close  to  the  maximum  amplitude  in  the  signal,  but  at  RKON  for  KASSERI 
the  maximum  amplitude  occurs  4 cycles  into  the  record.  (The  RKON  signals 
vary  substantially  in  character  as  a function  of  epicentral  location  and  many 
look  more  like  signals  at  HNME,  which  have  a simple  shape  for  all  events.) 
Though  these  maximum  amplitudes  have  been  measured  routinely  in  this  study, 
they  are  ignored  in  order  to  analyze  portions  of  the  signal  that  may  be  pre- 
dicted with  planar,  parallel-layer  earth  models:  these  4-  to  6-sec  late  arri- 
vals may  correspond  to  triplications  or  to  mantle-multipaths  that  cannot  pre- 
sently be  predicted.  Measurements  were  also  made  for  stations  WH2YK  and  FNWV, 
and  they  have  been  included  in  the  data  base.  However,  the  data  are  sparse 
because,  unlike  HNME  and  RKON,  the  stations  WH2YK  and  FNWV  were  out  of  opera- 
tion for  long  periods;  therefore,  measurements  from  those  2 stations  have  not 
been  analyzed. 

Figures  3a  through  3f  (from  Blandford,  1976)  show  the  variation  of  first 
motion  (labelled  "a"),  the  maximum  peak-to-trough  amplitude  (labelled  M£"), 
and  the  maximum  peak-to-trough  amplitude  divided  by  period  and  corrected  for 
system  response  at  that  period  (labelled  m£/T,  corrected  for  instrument  re- 
sponse") as  a function  of  yield,  for  several  values  of  the  absorption  para- 
meter t*  = T/Q,  where  T is  the  travel  time  and  Q is  the  average  quality  fac- 
tor along  the  path.  The  RDP  used  in  Figures  3a,b,c  is  for  granite  (von  Seg- 
gern  and  Blandford,  1972)  and  varies  as  a function  of  yield.  Although  this 
potential  is,  strictly  speaking,  for  granite  and  not  for  tuff,  which  is  more 
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Figure  3b.  Theoretical  amplitude-yield  curve  for  t*  * 0,  0.1,  0.2,  0.4,  O.n; 

granite,  amplitude  ol  thirst  motion,  with  surface  reflection,  pi' 
delay  equal  to  0.12Y  I sec,  with  Y in  kt . From  Blandford  ( 1^7(0 


10  * 10°  101  to2  to3  JO4 

YIELD  kt 

Figure  3d.  Theoretical  amplitude-yield  curves  for  t*  = 0,  0.1,  0.2,  0.4,  O.b; 

granite,  one-half  maximum  peak-to-trough  ampji^ude  of  signal,  with 
surface  reflection.  pP  delay  equal  to  0.12Y1/  second,  with  Y in 
kilotons.  From  Blandford  (197h). 
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Figure  Je.  Theoretical  amplitude-yield  curves  tor  t*  - 0,  0.1,  0.2,  0.-*,  0. 

granite,  one-half  maximum  peak-to-trough  amplitude  ot  signal, 
corrected  for  Instrument  response  at  measured  period,  T,  and 
divided  bv  T;  with  no  surface  reflection.  From  Blandford 


H 


A/T  (CORRECTED  FOR  PERIOD) 


YIELD  kt 

Figure  3f . Theoretical  amplitude-yield  curves  for  t*  = 0,  0.1,  0.2,  0.4,  0.6; 

granite,  one-half  maximum  peak-to-trough  amplitude  for  signal, 
corrected  for  instrument  response  at  measured  period,  T,  and 
divided  by  T,  with  surface  reflection.  pP  delay  equal  to  0.12Y1/3 
sec,  with  Y in  kt.  From  Blandford  (1976). 
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typical  of  NTS,  many  workers  such  as  von  Seggern  and  Klandtord  (1972),  Nopo- 
nen  (19/5),  Lyuke  et  al.  (19/6),  have  found  it  valid,  as  iar  as  spectral 
shape  is  conerned,  for  several  different  types  oi  competent,  saturated  rock 
observed  at  teleseismic  distances.  Blandford  (19/6)  concluded  that  although 
tuff  potential  derived  from  KATNTF.R  was  satisfactory  for  observations  at 
KNUT  and  MNNV  2 to  1 decrees  away,  it  was  inadequate  lor  explaining  teleseismic 
observations,  implying  that  the  RDP  was  a function  o!  take-oil  angle.  In 
this  study,  the  granite  potential  is  sufficient  to  explain  teleseismic 
observations  of  explosions  in  saturated  tuff. 

Figures  3a,  c,  and  e give  results  for  the  case  where  theoretically 

there  is  no  tree-surface  reflection  (pP).  This  could  occur,  for  example, 

if  the  free  surface  were  significantly  "rough"  or  if  there  were  significant 

spalling,  resulting  in  "trapping"  the  pi*  signal.  Figures  lb,  d , I give 

results  for  the  case  in  which,  in  theory,  there  is  perfect  tree-surface 

reflection.  In  this  case,  a theoretical  expression  must  be  found  lor  the 

reflection  delay  time,  given  hv  twice  the  depth  divided  by  the  average  velocity. 

Assuming  a velocity  ot  2.75  km/sec,  and  a containment  depth  formula  (V  in 
1/3  1/3 

kt)  of  h = 0.16Y  km,  yields  r = 0. 12Y  sec.  Many  events  at  Yucca  have 

1/3 

del. tvs  more  in  accordanct  with  0.2Y  , while  shots  in  granite,  such  as 

1/3 

PILEDRIVER  and  HARDHAT,  are  more  like  0.05Y  . Observed  data  should  lie 

between  the  curves  tor  ported  reflection  and  no  reflection.  Blandford 
(1976)  concluded  that  the  reflection  coefficient  for  pi*  at  NTS  for 
frequencies  greater  than  0.5  Hz  was  between  0.5  and  0.2.  Therefore,  in 
this  study,  the  experimental  amplitude  data  is  compared  to  Figures  3a,  c,  e, 
which  neglect  the  effect  of  pP. 
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ANALYSIS  OF  RKON  AND  HNME  DATA 


RKON 

Figure  4a  shows  a plot  of  the  base-10  logarithm  of  RKON  ji-phase  ampli- 
tude versus  Ms  from  von  Seggern  (1976).  For  the  smallest  events,  von  Seggern 
used  MNNV  to  determine  Mg.  For  larger  events,  and  for  events  after  the  1970 
closing  of  MNNV  as  an  LRSM  station,  von  Seggern  used  RKON.  Data  points  are 
missing  for  several  events  after  1970  because:  (1)  the  LR  was  too  weak  to  be 
measured  at  RKON;  (2)  because  the  RKON  short-period  data  was  unavailable;  or 
(3)  because  the  RKON  first-motion  amplitude  was  judged  to  be  severely  contam- 
inated by  noise. 

Figure  4a  also  includes  the  least-squares  line  with  fixed  slope  of  1.0. 
Figure  4b  shows  a plot  of  the  residuals  derived  from  Figure  4a.  This  is 
standard  procedure  in  this  paper,  except  when  plotting  period  T versus  Ms,  in 
which  case  both  the  slope  and  intercept  are  estimated.  Reference  to  Figure  1 
(use  transparent  overlay)  and  Table  I will  identify  which  event  is  plotted  at 
each  point  in  Figure  4b.  If  the  magnitude  residual  lies  between  -0.1  and 
+0.1,  then  a "0"  is  plotted  on  Figure  4b.  If  the  residual  is  0.1  to  0.2,  a 
"1"  is  plotted;  if  0.2  to  0.3,  a "2";  etc.  If  the  residual  is  -0.1  to  -0.2, 
then  an  "A"  is  plotted  on  Figure  4b.  If  the  residual  is  -0.2  to  -0.3,  a "B"; 
etc.  Note  that  in  Figure  4b,  the  amplitudes  from  events  in  the  western 
three-quar ters  of  Pahute  Mesa  are  substantially  smaller  than  the  amplitudes 
from  events  in  the  rest  of  NTS.  Amplitudes  from  events  in  the  eastern 
one-quarter  are  also  smaller  than  those  from  the  rest  of  NTS,  but  not  by  as 
great  a margin.  For  this  reason,  Pahute  Mesa  is  divided  into  two  regions: 
the  western  region,  for  whose  events  the  symbols  are  plotted  within  with 
an  open  circle  (0);  and  the  eastern  region,  where  events  are  plotted  within 
an  open  circle  with  a stick  at  "2  o'clock".  Events  shot  at  Rainier  Mesa 
are  plotted  within  a square  (□)  and  events  shot  at  Yucca  Flats  are  plotted 
within  a triangle  (A). 

In  Figure  4a,  there  are  also  plotted  two  curves  of  the  logarithm  of 
first-motion  amplitude  versus  yield;  these  are  from  Figure  3a,  with 
t*  * 0.4.  The  upper  line  represents  the  data  of  Rainier  Mesa  and  Yucca 
Flats  and  the  lower  line  represents  the  Western  Pahute  Mesa  data. 
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The  t*  value  Is  close  to  the  value  determined  by  Der  and  McElfresh  (1976)  for 
the  NTS-RKON  path. 

The  solid  line  trending  north-south  through  Yucca  Flats  in  Figure  4b  re- 
presents the  surface  trace  of  Yucca  Fault,  taken  from  detailed  geologic  maps 
(Johnson  and  Hibbard,  1957;  or  Fernald  et  al.,  1968).  Fernald's  map  shows 
that  the  hard  rock  at  Yucca  Flats  is  covered  by  a layer  of  superficial  de- 
posits varying  in  depth  from  700  m near  COMMODORE  (I)  to  0 m near  BOURBON 
(G). 

To  explain  the  large  variance  in  residuals  at  Yucca,  velocity  and  den- 
sity profiles  are  used  for  the  event  OBAR,  together  with  a Haskell  matrix  ap- 
proach. Results  show  that  for  an  event  just  below  the  surficial  deposits 
— paleozoic  interface,  the  teleseismic  amplitude  may  be  reduced  by  a factor 
of  2 or  more  by  cancelling  the  direct  wave  by  the  "internal  pP" 
wave  reflected  from  the  interface.  To  determine  if  this  type  of  calculation 
could  explain  the  variation  observed  at  Yucca  (shown  in  Figure  4b),  special 
calculations  can  be  performed  for  several  events  at  Yucca  and  then  compared 
with  the  predicted  and  the  observed  relative  amplitudes.  As  of  this  writing, 
these  calculations  have  not  been  performed. 

Missing  also  is  a convincing  explanation  for  the  low  amplitude  levels  in 
west  Pahute  Mesa.  However,  one  explanation  might  be  that  if  a large  high-Q 
plug  beneath  the  Silent  Canyon  Caldera  is  centered  on  event  GREELEY  (E)  (sug- 
gested by  Spence  (1973)  and  by  the  Caldera  surface  perimeter,  taken  from  the 
map  by  Orkild  et  al.  (1969)  and  drawn  in  Figure  1),  then  the  larger  first- 
motion  values  could  be  interpreted  for  STILTON  (w)  and  for  events  in  east 

Der,  Z.  A.,  and  T.  McElfresh  (1976).  The  effect  of  attenuation  of  the  spec- 
tra of  P waves  from  nuclear  explosions  in  North  America;  SDAC-TR-76-7, 
Teledyne  Geotech,  Alexandria,  Virginia.  DDC:  ADA  030857. 

Johnson,  M.  S.,  and  D.  E.  Hibbard  (1957).  Geology  of  the  Atomic  Energy  Com- 
mission Nevada  Proving  Grounds  Area,  Nevada;  U.S.  Geological  Survey  Bul- 
letin 1021-K,  U.S.  Government  Printing  Office,  Washington,  DC. 

Fernald,  A.  T.,  G.  S.  Corchary,  and  W.  P.  Williams  (1968).  Surficial  geolo- 
gic map  of  Yucca  Flats,  Nye,  and  Lincoln  counties;  U.S.  Geological  Sur- 
vey Map  1-550. 

Orkild,  P.  P.t  K.  A.  Sargent,  and  R.  P.  Snyder  (1969).  Geologic  map  of  Pah- 
ute Mesa  Nevada  Test  Site  and  vicinity,  Nye  County,  Nevada;  U.S.  Geolo- 
gical Survey  Map  1-567. 
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Pahute  as  showing  that  they  are  not  affected,  on  the  path  to  RKON,  by  the  ef-  I 

fects  of  the  high-Q  plug.  If,  however,  the  h 1 gli— Q plug  is  used  to  explain  1 

the  low  first-motion  amplitudes,  the  effect  must  be  due  to  ray  refraction  and  j 

divergence.  We  shall  see  that  the  variation  is  not  as  large  at  HNffF,  possibly  j 

because  refractions  at  the  edges  of  a high-Q  plug  do  not  serlouslv  affect  a 
more  vertically  departing  ray. 

Event  TYBO  (V)  shows  the  largest  deviation  from  the  average  M and  m 

s 

shown  in  Figure  4a.  In  a personal  communication,  L.  Ramspott  stated 
that  TYBO  had  anomalously  large  amplitudes  at  Las  Vegas  and  Mina  , Nevada,  to 
the  east,  and  anomalously  small  amplitudes  at  3 other  Lawrence  Livermore 
seismic  stations  at  other  azimuths.  In  Figure  4b,  the  contours  which  can  be 
drawn  for  west  Pahute  Mesa  suggest  that  if  similar  contours  were  drawn  for 
any  station  of  interest,  an  improved  amplitude  prediction  could  be  made  be- 
cause the  contours  show  TYBO  to  fit  the  overall  pattern  for  west  Pahute  Mesa. 

Note  that  because  these  contours  cannot  be  drawn  at  HNME,  this  smooth, 
dramatic  variation  in  contours  is  clearly  a source-station  effect  for  RKON, 
possibly  due  to  the  greater  angle  from  the  vertical  for  rays  departing  to 
RKON.  This  greater  angle  might  enhance  the  refraction  effects  as  the  rays 
leave  the  high-Q  plug. 

Plots  similar  to  those  in  Figures  4a, b have  been  created  with  M,,  re- 
placed with  log(yield).  Exactly  similar  conclusions  may  be  reached,  as  would 
be  expected  from  the  exceedingly  small  M :yield  scatter  found  by  von  Seggern 
(197b)  for  NTS  events  below  the  water  table.  The  low  variance  of  von  Seg- 
gern's  data  can  be  represented  by  a plot  of  Ms:yield  for  the  shots  with  de- 
classified yields,  given  in  Table  l.  These  data  are  well  represented  by  the 
equation  Ms  - log(yield)  + 2.4.  (von  Seggern  set  the  absolute  level  ot  Ms, 
using  a world-wide  single-station  network  (WWSSN)  of  10  stations  for  several 
large  NTS  shots.)  The  only  noticeable  change  is  that  events  above  the  water 
table  move  to  the  right  in  the  m^log  (yield)  plot  analogous  to  Figure  4a. 

This  shift  is  expected  because  these  events  do  not  couple  as  well  as  events 
below  the  water  table.  Figure  1 indicates  the  distance  (in  meters)  to  the 
water  table  for  events  above  the  water  table. 
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PILEDRIVER  ()),  a shot  in  granite,  has  a higher  amplitude  than  west  Pah- 
ute  events  possibly  because  of  the  source  medium.  Figure  4b  shows  that  the 
location  of  PILEDRIVER  also  can  be  used  to  explain  this  observation. 

Figure  5a  is  a plot  of  log(a)  versus  log(£)  at  RKON,  The  scatter  shows 
substantial  variation  in  wave  shape.  Variations  appear  during  comparison  of 
the  waveforms  for:  SCOTCH  (J),  COMMODORE  (1),  REX  (5),  and  TYBO  (v)  (Figure 
5d).  SCOTCH  and  REX  are  very  similar  in  wave  shape  and,  like  TYBO,  have  a 
small  first  motion.  The  small  first  motion  is  typical  of  Pahute  Mesa  events. 
Yucca  events,  however,  have  a relatively  large  first  motion  that  is  indicated 
by  the  short-period  vertical  trace  for  COMMODORE. 

Figures  5a, b show  a substantial  difference  in  the  log(a^)  versus  log(c) 
ratios  between  Pahute  and  Yucca.  This  effect  could  stem  from  multipathing 
or  from  the  fact  that,  because  the  average  velocity  in  the  medium  is  lower 
at  Yucca  than  at  Pahute,  the  pP  interference  with  the  c_  amplitude  is  greater 
at  Yucca,  so  that  the  first  motion  of  the  signal  appears  large  in  relation  to 
the  rest  of  the  signal. 

To  test  this  hypothesis,  the  log(a^)  residual  versus  the  theoretical  2- 
way  travel  time  (2x  depth/average  velocity)  is  plotted  (Figure  5c) . The  tra- 
vel time  has  been  normalized  by  the  dominant  period  which  the  received  P- 
wave  signal  would  have  if  not  affected  by  pP.  (This  dominant  period,  a func- 
tion of  yield,  is  determined  later  in  this  report.)  The  result  shown  in  Fig- 
ure 5c  indicates  a definite  trend,  that  suggests  that  pP  does  have  an  effect. 
Whether  the  trend  follows  predicted  detailed  theory  can  only  be  determined 
by  detailed  calculations  of  complete  theoretical  waveforms,  using  techniques 
discussed  by  Blandford  (1976). 

One  possible  physical  explanation  for  the  behavior  plotted  in  Figure  5c 
is  that  the  a amplitude  is  actually  unaffected  by  pP  in  any  of  the  measure- 
ments. However,  for  very  short,  scaled  delay  times,  the  c_  phase  is  substan- 
tially reduced  by  the  reflection  of  £.  Thus,  a has  a positive  residual  with 
respect  to  the  expected  value  for  £.  At  slightly  larger  delays,  the  £ phase 
is  enhanced,  while  for  even  larger  delays  it  is  again  reduced.  The  fact 
that  the  a minimum  (£  maximum)  occurs  at  a time  delay  of  approximately  0.5 
times  the  dominant  period  is  consistent  with  this  interpretation. 
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The  empirical  results  In  Figure  Sc  surest  that  the  c amplitude,  which 
in  most  cases  will  he  the  phase  measured  tor  magnitude,  may  varv  by  at  least 
0.r>  units  of  magnitude  due  to  depth  ot  burial.  While  amplitude  can  In- 
crease by  only  0.1  magnitude  units,  it  may  decrease,  theoretically,  to  zero. 
Thus,  the  range  indicated  in  Figure  Sc  does  not  seem  implausible.  In  this 
report,  the  minimum  delay  occurred  tor  shots  in  granite;  tor  a velocity  ot 
s.  S km/sec,  the  delay  for  Pll.KOKIVKR  is  0.17  sec. 

Figures  6a, b show  plots  ot  logic)  versus  Ms  anil  the  residuals  ot  log(_c) 
witli  respect  to  latitude  and  longitude.  In  this  case,  variations  attributed 
to  multipathing  or  to  ray  convergence  and  divergence  in  Figures  7a, b would  be 
confounded  by  the  variations  due  to  the  pi’  Interference  discussed  in  connec- 
tion with  Figures  5a, b.  Thus,  simple  contours,  like  In  Figure  4b,  cannot  hi1 
drawn  in  west  I’ahute.  Curvature  as  a function  ot  yield  tor  the  mean  line  for 
west  i’ahute  seems  more  ot  a possibility  in  Figure  6a  than  in  Figure  4a.  Indi- 
cations strongly  suggest  that  it  the  waves  ot  KKON  have  a short  path  (or  no 
path)  through  the  root  of  the  Silent  Canyon  Caldera,  then  the  c amplitude  is 
substantially  Increased.  However,  this  effect  may  be  confounded  by  the  mag- 
nitude-yield curvature.  For  example,  though  event  STILTON  (w)  has  a clear 
path  from  its  position  northwest  of  the  caldera  to  the  station  at  KKON,  its 
yield  is  low  enough  to  make  its  corner  frequency  high,  which  will  result  in 
a relatively  large  magnitude.  Note  in  Figure  6b  that  SCOTCH  (.1)  and  CAMFM- 
RFRT  (z)  appear  to  be  in  a transition  region  between  east  and  west  Pahute. 

Figure  7 shows  a plot  of  log(£>,  corrected  for  pP,  versus  log(a).  The 
correction  factor  is  derived  from  the  line  through  the  data  in  Figure  Sc,  and 
is  done  to  obtain  a smoother  m^:Ms  plot  by  removing  the  perturbing  effects,  of 
pP . The  result  should  be  similar  to  Figure  4a  — log (a)  versus  M — and 
this  indeed  seems  to  be  the  general  result. 

In  Figure  8a  the  period  of  the  c phase  at  KKON  is  plotted  as  a function 
ot  c versus  . ['he  plot  shows  a strong  trend  over  the  depicted  range  ot 
period  (about  1.0  to  1.5  sec)  and  range  of  M (about  2.7  to  5.5)  — roughly 
equal  to  a yield  ranging  from  less  than  5 kt  to  slightly  more  than  l megaton. 
This  variation  of  period  with  M was  used  in  Figure  5c  to  determine  the  dom- 
inant period  with  which  to  normalize  the  theoretical  2-way  pP  delay  time. 

The  same  relationship  will  be  used  in  similar  plots,  such  as  Figure  8c,  below. 
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Figure  6b.  Residuals  of  log  (c)  with  respect  to  the  dashed  line  in  Figure 
6a  as  a function  of  latitude  and  longitude  at  NTS.  Events  can 
be  identified  by  comparison  with  Figure  1 and  Table  I. 
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Note  the  large  scatter  about  the  mean  Line  In  Figure  Ha,  wliieli  the  au- 
thors contend  stems  from  geographic  location  and  depth  of  burial  (via  the  ef- 
fect of  pP  Interference).  Figure  Fib  shows  the  residuals  ibout  the  regression 
line  plotted  as  a function  of  location;  the  contrast  between  east  and  west 
Pahute  is  apparent.  Again,  events  SCOTCH  (.1)  and  CAMKMHKKT  (/.)  seem  to  be  in 
a transition  region.  In  Figure  He  these  period  residuals  are  plotted  as  a 
function  of  of  normalized  2-wav  pP  delay  time,  and  a clear  trend  is  evident. 

The  events  at  either  end  of  the  trend  are  PILKDKIVER  ()),  and  COBBLER  (P). 
PILKDKIVER  has  an  unusually  short  scaled  delay  time  because  ol  the  high  com- 
pression wave  velocity  in  granite,  and  COBBLER  has  a longer  than  expected 
delay  time  because  it  was  overburled  for  its  yield  and  the  medium  had  a low 
compression  wave  velocity.  This  graph  suggests  that  the  effects  of  pP  alone 
may  change  the  dominant  period  by  as  much  as  0.7  sec.  Note  that  this  period 
correction  will  have  an  effect  on  m^.  At  constant  amplitude  on  a short-per- 
period  LRSM  system,  this  possible  *0.  15-sec  variation  in  period  around  1.2  sec. 
corresponds  to  a range  of  0.14  magnitude  units  in  the  factor  1/CT.  In  sum- 
mary, Figure  8a  shows  the  effect  of  yield  on  the  dominant  period.  Figure  8b 
shows  the  effect  of  geographic  location,  and  Figure  He  shows  the  effect  of  pP. 

While  all  effects  are  significant,  note  that  they  are  linked;  for  example, 
large  events  are  deep  and  in  Pahute  Mesa.  Only  calculations  of  each  event 
can  disentangle  these  effects. 

Most  of  the  trend  toward  short  periods,  shown  in  Figure  He,  can  be 
traced  to  the  2 granite  shots.  Although  the  periods  mav  be  short  because 
of  the  unique  material  properties  of  granite,  that  argument  is  weak  because 
the  granite  RDP  source  seems  to  be  adequate  to  explain  the  teleseismic 
magnitude  versus  yield  curves  for  well-contained  explosions  in  all  source 
regions  at  NTS.  Thus,  there  is  no  clear  reason  to  omit  the  data  from  the 
shots  truly  set  off  in  granite. 

In  Figures  9a,  h,  are  plots  similar  to  those  in  Figures  ba,h.  The  differ- 
ence is  that  in  Figures  9a, h the.  £ amplitude  has  been  corrected  for  system 
response  at  its  observed  period,  and  the  resulting  amplitude,  divided  by  T, 
has  been  plotted.  Considering  the  variety  of  Influences  >n  observed  period 
already  discussed,  probably  no  correction  for  period  woul 1 have  any  substan- 
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tial  clear-cut  effect  on  the  scatter  of  the  mb:Ms  diagram;  in  fact.  Figures 
6 and  9 are  remarkably  alike.  However,  note  that  both  theoretically  (without 
pP)  and  experimentally,  log(c/T)  has  a slope  closer  to  1.0  near  1000  kt  than 
does  log (c) . 

HNME 

Figures  like  the  series  just  studied  from  RKON,  but  with  data  from 
HNME  will  now  be  examined.  From  consideration  of  values  derived  in  Der  and 
McElfresh  (1976),  a value  of  t*  - 0.6  is  used.  In  Figures  10a, b like 
Figures  4a. b there  is  evidence  of  amplitude  differences  between  source  re- 
gions at  NTS.  Because  the  signal-to-noise  ratio,  especially  for  first  motion, 
is  lower  at  HNME,  the  scatter  for  low  magnitudes  (Figure  10a)  is  probably 
greater  than  it  would  otherwise  be.  The  events  with  question  marks  after 
their  symbols  (Figure  10b)  have  especially  questionable  first-motion  ampli- 
tudes due  to  noise. 

Figure  11a  shows  a plot  of  log  (a)  versus  log(jc)  at  H'JME.  Compared  to 
the  equivalent  plot  for  RKON  (Figure  5a),  the  HNME  waveform  is  only  weakly 
affected  by  location  within  NTS.  With  respe«  t to  amplitude,  however,  by 
plotting  log(c)  versus  Mg  (Figures  12a, b),  it  again  appears  possible  to  dis- 
tinguish between  west  Pahute  and  the  remainder  of  NTS. 

The  same  patterns  of  results  are  obtained  for  the  c-phase  amplitudes 
at  both  RKON  and  HNME,  a regional  and  a teleseismic  station,  suggesting 
that  network  averages  would  show  the  same  effects.  Note,  however,  that 
RKON  and  HNME  are  similar  in  distance  and  azimuth,  and  that  the  average 
results  from  a well-distributed  world-wide  network  might  reduce  the  above- 
mentioned  effects  if  they  were  caused  by  multipathing  in  the  upper  crust 
structure. 

Note  in  both  Figure  10a  and  12a,  the  west  Pahute  data,  signified  by 
open  circles,  forms  a curved  line  that  is  much  better  defined  than  the 

line  for  the  entire  population. 


function  of  latitude  and  longitude  at  NTS.  Events  may  be  identi- 
fied by  comparison  with  Figure  1 and  Table  I. 
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tude  at  NTS.  Events  may  be  identified  by  comparison  with 
Figure  1 and  Table  I. 


Figure  11a.  a as  a function  of  c at  HNME.  The  dashed  line  is  the  least- 
squares  line  of  slope  1.0  through  all  the  data. 
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Figure  12b-  .c  residuals  with  respect  to  the  dashed  line  in  Figure  12a  at 
HNME  as  a function  of  latitude  and  longitude  at  NTS.  Events 
may  be  identified  by  comparison  with  Figure  1 and  Table  I. 


Using  tlu*  sol  ill  l i ih*  in  Killin'  I U*  to  correct  tin*  c phase  results  in 
Figure  11,  which  is  not  significantly  il  liferent  I rum  12a.  This  result  i m 
not  surprising  considering  the  small  it).  1 magnitude  units)  maxlmum-to- 
minimum  excursion  curve  drawn  in  Figure  lie. 

In  Figure  14a,  like  Figure  8a,  is  a variation  in  dominant  period  between 
1.0  and  1 . t>  seconds  as  yield  varies  between  5 and  1000  kilotons.  The 
percentage  variation  resulting  I rom  oplcontral  location  and  depth  (Figures 
I4b,c)  seems  somewhat  less  for  HNMF.  than  lor  KKON.  The  regression  line 
drawn  in  Figure  14c  (the  same  drawn  in  Figure  8c)  is  perhaps  too  steep.  On 
tin*  other  hand,  unlike  KKON  the  trend  in  Figure  8c  does  not  seem  dependent  on 
the  granite  point. 

Inspection  ot  Figures  8a  and  14a  leaves  little  doubt  that  a variation 
exists  in  the  relation  ol  pel iod  to  vield.  Springe!  and  Hannon  (1021)  con- 
cluded that  there  was  no  variation  at  teleseismic  distances.  However,  they 
had  available  tar  fewer  events  to  plot  and  were  torced  to  accept  the  period 
measurements  made  in  routine  shot  analyses.  For  the  larger  amount  ol  data 
analyzed  in  this  report,  the  scatter  is  so  great  th.it  select  Ion  ol  only  It) 
points  (tlii*  number  ol  data  available  to  Springer  and  Hannon)  could  easily 
lead  to  a conclusion  that  there  was  no  trend  to  the  data.  Also,  many  ot  the 
routine  shot  measurements  are  undoubtedly  contaminated  by  "talse  cycles",  a;, 
discussed  in  connection  with  Figure  2. 

Figure  15a  shows  that  correcting  for  amplitude  response  at  the  measured 

period  T and  dividing  by  T results  in  little  change  in  the  scatter  ol  the  m : 

M curve.  However,  note  that  compared  with  Figure  12a,  the  experimental 
s 

log  (c/T)  curve  as  a slope  closer  to  1.0  than  does  the  log  (c)  curve. 
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gure  15b.  HNME  Station:  Residuals  of  log  (c/T)  with  respect  to  the 

dashed  line  in  Figure  15a  as  a function  of  latitude  and  longi- 
tude at  NTS.  Events  are  identified  in  Table  I and  Figure  1. 
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SUMMARY  AND  DISCUSSION 

In  this  report,  von  Seggern's  (1976)  near-equivalence  of  log(yield)  and 
Ms  has  been  used  to  show  that  both  period  and  amplitude  of  short-period  tele- 
seismic  body  waves  vary  with  yield,  effects  of  pP,  and  geographic  location. 

The  variation  from  direct  proportionality  of  amplitude  versus  log(yield)  at 
megaton  levels  is  on  the  order  of  0.1  to  0.2  magnitude  units,  and  the  effects 
of  pP  have  been  shown  to  result  in  a spread  of  about  0.34  magnitude  units. 

The  variation  with  geographic  region  at  an  individual  station  is  different 
for  RKON  and  HNME,  but  a median  value  comparing  west  Pahute  Mesa  with  Yucca 
Flats  might  be  taker  as  0.4  magnitude  units. 

The  period  varies  from  approximately  1.0  sec  for  a 5 kt  explosion  to  1.6 
sec.  for  a 1000  kiloton.  The  period  also  shows  a substantial  variation  due 
to  pP  effects  and  local  geology.  With  these  facts  in  mind,  one  would  not 
expect  it  would  be  advisable  to  correct  measured  amplitudes  for  the  system 
response  at  the  measured  period  if  all  measurements  were  made  on  a common 
system.  The  problem  is  made  worse,  of  course,  if  one  considers  that  period 
measurements  are  notoriously  difficult  and  that  the  signal  is  not  monochroma- 
tic. 

As  reported  in  this  study,  magnitude  and  period  effects  of  the  size  ob- 
served at  single  stations  suggest  that  routine  log(c/T)  measurements  cannot  be 
adequate  to  the  task  of  estimating  seismic  magnitude  if  a testing  organiza- 
tion varies  the  depth  of  the  events  to  an  appreciable  degree.  Since  the 
depth  effect  would  be  common  to  all  stations  in  a network,  the  resultinv  Mac 
could  not  be  reduced  by  network  averaging. 
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